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SUMMARY 


The  effects  of  transverse  shear  deformations  upon  the  classical  buckling 
load  of  fiber-reinforced  cylindrical  shells  under  uniform  axial  compres¬ 
sion  have  been  analyzed.  A  stability  determinant  which  upon  evaluation 
yields  an  expression  for  the  buckling  load  has  been  developed  using  a 
modified  form  of  the  Reissner  variational  principle.  The  buckling  loads 
predicted  by  the  stability  determinant,  with  transverse  shear  effects 
neglected,  have  been  calculated  and  shown  to  agree  with  previously- 
published  results  for  boron-epoxy  and  glass-epoxy  cylinders.  Inclusion 
of  the  transverse  shear  effects  in  the  two  cases  investigated  shows 
little  reduction  in  the  classical  buckling  loads.  For  general  appli¬ 
cation,  charts  are  presented  to  give  stability  criteria  for  fiber- 
reinforced  composites  as  a  function  of  the  geometric  and  mechanical 
properties.  Effective  transverse  shear  rigidities,  established  by  experi¬ 
ment,  must  be  developed  in  order  to  estimate  realistically  the  reduction 
in  buckling  load  from  the  charts  presented.  For  the  determination  of 
natural  frequencies  of  fiber-reinforced  shells,  with  transverse  shear 
effects  included,  a  frequency  equation  in  determinant  form  is  presented. 
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FOREWORD 


The  work  reported  herein  constitutes  a  portion  of  a  continuing  effort  be¬ 
ing  undertaken  at  Stanford  University  for  the  U.  S.  Army  Aviation  Materiel 
Laboratories  under  Contract  DAAJ02 -68 -C-003 5  (Task  1F16220UA17002)  to 
establish  accurate  theoretical  prediction  capability  for  the  static  and 
dynamic  behavior  of  aircraft  structural  components  utilizing  both  con¬ 
ventional  and  unconventional  materials.  Previous  contracts  supported 
investigations  which  led,  in  part,  to  the  results  presented  in  references 
2,  3,  k,  and  5. 
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INTRODUCTION 


Due  to  their  high  strength- to-weight  ratios,  fiber-reinforced  cooposites 
are  increasir-ly  be^Pg  used  for  structural  applications.  In  the  anal¬ 
yses  to  determine  the  structural  characteristics  of  fiber-reinforced 
beams,  plates,  and  shells,  the  structural  members  have  been  assumed  to 
be  infinitely  rigid  with  respect  to  transverse  shear  deformations.  The 
objective  of  the  present  analysis  is  to  investigate  the  effect  of  re¬ 
laxing  the  rigid-transverse-shear  assumption. 

The  usual  fiber- reinforced  composite  structure  is  laminated  from  thin 
layers,  each  having  unidirectional  fiber  reinforcement.  For  a  given 
layer,  the  ratio  of  the  inplane  shear  rigidity  to  the  extensional 
rigidity  is  much  lower  than  the  corresponding  ratio  in  a  homogeneous 
isotropic  material.  Therefore,  transverse  shear  deformations  during 
bending  in  laminated  fiber-reinforced  structures  can  be  expected  to  be 
much  more  significant  than  shear  deformations  in  homogeneous  isotropic 
structures.  In  the  present  study,  an  analysis  which  includes  the  effect 
of  transverse  shear  deformation  on  the  elastic  stability  of  fiber- 
reinforced  circular  cylindrical  shells  subjected  to  uniform  axial  com¬ 
pression  is  developed.  Applications  are  made  to  the  special  cases  of 
boron-epoxy  and  glass-epoxy  materials. 

The  analysis  is  based  upon  the  use  of  the  Reissner1  variational  prin¬ 
ciple.  First,  the  indirect  variational  procedure  is  used  to  establish 
the  governing  equations  and  the  general  boundary  conditions  for  the 
problem.  Then,  to  calculate  buckling  loads,  the  direct  variational 
approach  (Rayleigh-Ritz  method)  is  applied.  The  general  buckling  cri¬ 
teria  results  are  given  in  the  form  of  a  stability  determinant,  and  the 
specific  problem  solution  results  are  given  in  the  form  of  charts.  The 
nature  of  the  development  is  such  that  with  very  little  additional  effort 
the  natural  frequency  equation  for  fiber-reinforced  shells  with  trans¬ 
verse  shear  effects  included  can  be  obtained;  the  frequency  equation  is 
presented  in  an  appendix. 
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GENERAL  THEORY 


PROBLEM  STATEMENT 

The  general  problem  considered  here  is  the  investigation  of  the  effect 
of  transverse  shear  deformations  on  the  buckling  loads  of  circular 
cylindrical  shells  made  of  fiber-reinforced  composite  materials 
under  uniform  axial  compression.  The  .shell  geometry  and  coordinate 
system  used  are  defined  in  Figure  1.  The  fiber  reinforcement  within 
each  lamina  is  unidirectional;  however,  the  orientation  of  a  given  layer 
with  respect  to  adjacent  layers  can  be  arbitrary.  Also,  the  thicknesses 
of  each  of  the  layers  may  differ.  The  shell  wall  is  shown  in  Figure  2. 

REISSNER  FUNCTIONAL 

In  the  present  analysis,  a  variational  procedure  is  applied  to  the 

1  2  H  ; 

Reissner  form  of  the  potential  energy.  As  shown  by  Mayers  et  al. ’  , 
the  Reissner  formulation  and  modified  forms  thereof  are  demonstrated  to  be 
better  alternatives  to  the  potential  energy  principle  for  the  nonlinear 
analysis  of  plates  and  shells.  To  facilitate  the  extension  of  the  pres¬ 
ent  analysis  to  nonlinear  problems  and  to  preserve  continuity,  a  modi¬ 
fied  Reissner  principle  is  used  for  the  present  linear  analysis.  For  a 
thin  plate  or  shell  in  which  normal  stresses  in  the  direction  perpen¬ 
dicular  to  the  surface  are  neglected,  the  Reissner  functional  becomes 
(see  reference  l) 


U"  =  j  e  +  o-  e  +  r  .  +  t  -v  +  t  v  -  F '  |  dV 

./  I  x  x  y  y  xy  xy  xz  xz  yz  yz  1 

V 


(1) 


where  F'  is  the  complementary  energy  density  and  is  defined  by 


X  y  XZ  yz 

=  I  e  la  +  f  ev^°V  +  /  7  dr  +  f  y  dr  +  I  y  dr  (2) 

J  x  x  J  y  y  )  'xy  xy  /  xz  xz  /  yz  yz 
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xy 

o  o 


In  the  linear  elastic  case, 


F'  =  r  a  e  +  a  e  +  t  y  +  t  y  +t> 

2  x  x  y  y  xy'xy  xz'xz  yz7yz 


(3) 


2 


REFERENCE  SURFACE 
OF  THE  SHELL 


Figure  1.  Sign  Convention  and  Geometry. 


Figure  2.  Element  of  Shell  Wall. 
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RE  ISSUER  FUNCTIONAL  -  PLATES  AND  SHELLS 


The  strains  can  be  expressed  in  terms  of  the  reference-surface  strains 
and  the  relative  rotations  (corresponding  to  the  Kirchoff-Love  assump¬ 
tions)  as 
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(4) 


After  integration  through  the  thickness,  the  Reissner  functional  becomes 
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where  the  following  usual  definitions  of  the  stress  resultants  apply: 
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and  where  y  and  y  represent  average  transverse  shear  strains. 
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CONSTITUTIVE  LAW 

Discussions  of  the  linear-elastic  lav  for  laminated  anisotropic  plates 
and  shells  can  he  found  in  references  6,  7,  8  and  9.  Assumption  of  a 
state  of  generalized  plane  stress  gives  the  stress-strain  relations  for 
the  kth  layer  as 


where  the  symmetrical  matrix  [C'  '3  represents  the  elastic  coefficients 
associated  with  the  k  n  lamina.  The  elastic  coefficients  depend  upon 
the  elastic  properties  of  the  lamina  referred  to  a  set  of  elastic  axes 
and  the  orientation  of  the  elastic  axes  vith  respect  to  the  coordinate 
system  of  the  reference  surface  in  the  structure. 

The  stress  resultants  can  now  be  related  to  the  reference-surface  strains 
and  the  relative  rotations  and  twist  by 


where  .Aj  ,  LDJ  ,  and  lD  j  are  3x3  matrices  whose  elements 
depend  upon  the  physical  and  geometric  properties  of  the  laminae  (see 
Appendix  VI ).  Coupling  of  the  reference-surface  strains  and  the  loads 

with  the  bending  moments  and  the  relative  rotations  and  twist  due  to  the 
nonsymmetrical  arrangement  of  the  laminae  relative  to  the  reference  sur¬ 
face  is  produced  by  [D).  The  details  of  the  evaluation  of  the  ele¬ 
ments  of  the  various  matrices  are  given  in  references  6  through  9* 

The  transverse  shear  forces  and  average  transverse  shear  strains  are 
assumed  to  be  related  by 


3* 


hG  ^ 

X  xz 


and 


a  =  hG  7 
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(10) 


where  G  and  G  are  the  effective  transverse  shear  rigidities. 

x  y  c 

STRAIN-PIS  PLACEMENT  RELATIONS 

For  a  cylindrical  shell,  using  the  Donnell10  approximation,  the  reference 
surface  strains  for  small  deflections  are  given  by 
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In  a  manner  analogous  to  that  of  Libove  and  Batdorf11,  the  relative  ro¬ 
tations  and  the  relative  twist  are  given  by  the  relations 
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(12) 


in  which  the  terms  involving  derivatives  of  ~xz  and  7yZ  (the  average 
transverse  shear  strains)  are  corrections  to  the  usual  curvature  terms 
introduced  to  account  for  the  effects  of  transverse  shear  deformations. 


MODIFIED-REISSNER  FUNCTIONAL 

In  a  manner  similar  to  that  of  references  2,  3,  k,  and  5,  the  Reissner 
functional  of  equation  (5)  is  now  modified  so  that  the  force  resultants, 
the  reference  surface  displacements,  and  the  average  transverse  shear 
strains  are  the  only  variationally  independent  quantities.  The  bending 
and  twisting  moment  resultants  and  the  transverse  shear  force  resultants 
are  eliminated  from  the  Reissner  functional  through  the  use  of  appropri¬ 
ate  constitutive  relations. 


To  achieve  the  modification,  equation  (8)  is  inverted  and  substituted 
into  equation  (9).  with  the  result  that 
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where  [a]  ,  [d']  ,  [d]  ,  and  [d  ]  are  3x3  matrices  which  are  de¬ 
fined  in  terms  of  [A]  ,  [D]  ,  and  [D  ]  in  Appendix  VI.  The  comple¬ 
mentary  energy  per  unxt  surface  area  is  calculated  by  integrating 
equation  (3)  across  the  thickness;  by  applying  the  definitions  of 
equations  (6);  and  by  eliminating  the  transverse  shear  resultants, 
reference-surface  strains,  and  bending  moments  through  the  use  of 
equations  (10),  (13),  and  (14).  Such  a  procedure  -  the  same  as  that 


used  by  Khot~  for  evaluating  the  strain  energy  —  yields 
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The  potential  associated  with  the  applied  axial  compression  load  is 
given  by  (see  reference  12) 


(15) 
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Combination  of  equations  (5),  (10),  (ll),  (12),  (14),  (15),  and  (16) 
gives  the  following  modified  form  of  the  Reissner  total  functional  to  be 
used  in  the  variational  procedures  of  tha  remainder  of  this  analysis. 
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yz 
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ay2  ay  / 
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^  \  axay  ay  ax  / 


-  (dr,  N  +  d^_N  +  drr. 
'  61  x  62  y  66 


-  -  (a^N2  +  a.-N2  +  a^N2  +  2a,„N  N  +  2a. N  +  2a^B  N  ) 
„  11  x  22  y  66  xy  12  x  y  T.6  x  xy  26  y  xy 

2 


1 

+  — 
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xy  2b  y  xy 
2 


+4(^.5E)+Weiz..S-.3B) 

\  ax  ax  /  \  ay  ay  •  \  axay  ay  ax  / 
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(Continued) 


An  indirect  variational  approach  to  the  problem  result*  in  the  establish¬ 
ment  of  the  governing  equations  (Euler  equations)  end  the  associated 
boundary  conditions  of  the  surface.  The  details  of  the  indirect  approach 
are  given  in  Appendix  I;  the  resulting  Euler  equations  are 

Bn  bn 

-s  +  -s  „  o 

Bx  By 


(IB) 


9 


Dm  Dm 
_x  +  —EL 

Dx  Dy 


Dm  Dm 
—L  +  —EL 

Dy  ox 


(20) 


D2M  D2M  c2M 

+  2  _ 2SL  +  -JL 

Dx“  DxDy  Dy 


N 

.  JL 

R 


D2' 


w 


+  N 


x  s  2 
o  cx 


(21) 


which  apply  to  flat  (R  -4  <-o )  and  circularly  curved  plates  as  well  as 
to  circular  cylinders.  Equations  (18)  represent  equilibrium  in  the 
x-  and  y-directions,  respectively.  Equation  (19)  represents  stress- 
displacement  compatibility.  Equations  (20)  represent  moment  equilibrium 
about  the  x-  and  y-axes,  respectively.  Finally,  equation  (21)  repre¬ 
sents  equilibrium  in  the  direction  normal  to  the  reference  surface. 
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METHOD  OF  SOLUTION 


ASSUMED  SOLUTION  FUNCTIONS 

The  generalized  displacements  and  forces  satisfying  the  compatibility 
conditions  derived  in  Appendix  II  and  Appendix  III  and,  at  least,  the 
geometric  boundary  conditions  (simple  support  for  the  present  problem) 
are  assumed  to  be 


w  =  W  sin  sin 

L  b 


_  mnx  .  nny 

T  cos  —r—  sin  — t” 
x  L  b 


„  .  mnx  nny 

ry  Sln  T  C0S  b 


ux  sin  cos  211  +  u2  cos  52*  sin  22£ 


„  mnx _ nny  ^  ,  mnx  .  nny 

=  v1  sin  cos  +  v^  cos  -j-  sin 


■m  »r  ,  »,  •  mnx  .  nny 

N  =  N  +  N„  sin  — sin  -r*- 

x  x  3  L  b  ■ 


,T  «  •  mnx  •  nny 

N  =  N,  sin  —t—  sin  -r*- 
y  i  l  b 


„  ,T  mnx  nny 

Nxy  =  N2  C0S  IT  C0S  b 


where  w  *  Tx  ,  T  ,  ,  Ug  ,  vx  ,  vg  ,  ^  ,  N2  ,  and  Nj  become  the  in¬ 

dependent  quantities  in  the  direct  variational  procedure. 


STABILITY  DETERMINANT  DEVELOPMENT 

Substitution  of  the  assumed  functions  into  equation  (17 ),  integration 

over  the  prescribed  limits,  and  subsequent  variation  with  respect  to 

the  9  independent  quantities  yields  a  set  of  9  linear  homogeneous 

equations.  The  determinant  of  the  coefficients  of  the  homogeneous 

equations  establishes  the  condition  for  equilibrium  in  the  buckled 

state.  An  expression  for  the  applied  axial  load,  N  ,  required  for 

o 

equilibrium  in  the  buckled  shape  as  a  function  of  the  shell  geometric 

and  elastic  properties,  the  number  of  circumferential  waves  (N)  , 

and  the  buckle  aspect  ratio  (u)  results  when  the  determinant  is 

evaluated.  For  a  given  shell,  the  buckling  load  corresponds  to  the 

lowest  value  of  Nv  for  all  the  possible  combinations  of  N  and  u. 
xo 

The  stability  determinant  as  well  as  the  details  of  the  buckling  load 
calculation  Is  given  in  Appendix  IV. 
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RESULTS  AND  DISCUSSION 


The  primary  result  of  the  present  analysis  is  the  establishment  of  the 
stability  determinant  given  by  equation  (40)  in  Appendix  IV.  When  the 
shell  geometry  and  composite  material  construction  are  defined,  the 
axial  buckling  load  for  a  circular  cylindrical  shell  made  from  any 
number  of  arbitrarily  orientated  fiber-reinforced  layers  can  be  deter¬ 
mined  through  numerical  evaluation  of  the  stability  determinant. 

Specifically,  the  effects  of  transverse  shear  deformation  on  the  axial 

buckling  load  of  fiber-reinforced  shells  made  of  two  different  matrix- 

fiber  combinations  have  been  calculated  using  equation  (40).  Both  of 

the  shells  considered  have  a  6.0  in.  mean  radius  and  are  laminated  from 

0.012  in.  thick  layers  with  all  fiber  reinforcement  oriented  in  the 

longitudinal  direction.  The  lamina  elastic  properties  for  the  two 

ft  Q 

cases  studied  are  the  same  as  those  used  by  Khot  and  are  listed 
below. 


Boron-Epoxy  Composite 


Eu 

=  40.0  X  106 

lb/in.2  =  C1;L(1- 

V12V21> 

E22 

=  4.5  x  106 

lb/in.2  =  C22(l- 

V12V21> 

V12 

=  0.25 

v21  -  0.028 

G 

=  1.5  x  106 

lb/in. 2  =Ca 

Glass-Epoxy  Composite 


'll 

=  7.5  X  106 

lb/in.2  =  cu(1-v12v2l> 

'22 

=  3.5  x  106 

lb/in.2  =  c22(1-,i2v2i) 

12 

=  0.25 

V21  =  0.18 

G 

=  1.25  X  106 

ib/m.2  =  ca 

The  axial  buckling  loads  for  the  cases  of  infinite  shear  rigidity  (no 
transverse  shear  effect),  along  with  the  reductions  in  the  axial  buckling 
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loads  due  to  finite  transverse  shear  deformations,  are  shown  in 
Figures  3  and  4  for  boron-epoxy  and  glass-epoxy  composites,  respec¬ 
tively.  The  axial  buckling  loads  given  in  Figures  3  and  4  for  the 
cases  of  i-layer  shells  are  the  same  as  those  calculated  by  Khot  on 
the  basis  of  infinite  transverse  shear  rigidity  of  the  medium. 


REDUCTION  IN  BUCKLING  LOAD  (% 


NUMBER  OF  LAYERS 


Figure  3.  Stability  Criteria,  Including  Transverse  Shear  Effects,  for 
Boron-Epoxy  Cylindrical  Shells  (R  =  6.0  in.,  h.  =  0.012  in.) 


REDUCTION  IN 

BUCKLING  LOAD  (%)  BUCKLING  LOAD  (LB/IN.*  10 


NUMBER  OF  LAYERS 


Figure  4.  Stability  Criteria,  Including  Transverse  Shear  Effects,  for 
Glass-Epoxy  Cylindrical  Shells  (R  =  6.0  in.,  h  *=  0.012  in.) 
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CONCLUDING  REMARKS 


For  a  first  approximation  to  the  effects  of  transverse  shear  deformation 

on  the  buckling  load  of  axially  compressed  boron-epoxy  and  glass-epoxy 

fiber- reinforced  shells,  the  effective  shear  rigidities  (5  and  G  ) 

x  y 

can  be  assumed  to  have  values  between  the  shear  modulus  of  the  epoxy 

2 

matrix  (about  200,000  lb/in.  )  and  the  inplane  shear  modulus  of  the 

6  2  6 

composite  layer  (I.25  x  10  lb/in.  for  glass-epoxy  and  1.  *5  x  10° 

p 

lb/in.  for  boron-epoxy).  Figures  3  and  U  show  that  in  this  range  of 
values,  the  buckling  loads  calculated  with  transverse  shear  deformations 
included  are  less  than  7  percent  below  those  calculated  with  the 
assumption  of  no  transverse  shear  deformation.  Actual  effective  trans¬ 
verse  shear  rigidities  can  best  be  determined  by  experiment. 

Due  to  the  similarity  in  the  formulation  of  the  free  lateral  vibration 
problem  and  the  buckling  problem  to  which  most  of  the  preceding  pages 
have  been  devoted,  a  frequency  equation  which  may  be  used  to  establish 
the  effects  of  transverse  shear  deformation  on  the  natural  frequencies 
of  circular  cylindrical  fiber-reinforced  shells  has  been  developed  in 
Appendix  V.  For  given  materials,  natural  frequencies  can  be  established 
by  finding  the  eigenvalues  of  the  determinant  (equation  Uj)  in  a  pro¬ 
cedure  similar  to  that  used  for  the  buckling  problems. 
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APIENDIX  I 


EULER  EQUATIONS  AND  BOUNDARY  CONDITIONS 
DERIVED  FRCM  THE  VARIATIONAL  PRINCIPLE 

The  vanishing  of  the  first  variation  of  equation  (17 )  with  respect  to 

N  ,N  ,N  ,  u  ,  v  ,  w  ,  >  ,  and  y  gives 
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-  *ll\-Nx  ‘  Vy^y  '  ‘  W',  * 


(Continued) 
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OW  o?w 

-  Nx - 

o  :,x 


dxdy 


=  0 


(??) 


After  integration  by  parts,  regrouping  of  terms,  and  application  of 
equations  (ll)  and  (12)  along  with  the  constitutive  relations  (10) 
and  (14),  equation  (25)  becomes 
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(26) 


For  the  above  expression  to  vanish  for  the  specified  arbitrary  states 
of  stress  and  displacement  consistent  with  the  geometric  boundary  con¬ 
ditions,  each  of  the  terms  must  vanish  independently.  The  first  eight 
terms  of  equation  (26)  lead  to  the  Euler  equations  and  the  remaining 
terms  establish  the  boundary  conditions,  as  follows : 


Inplane  equilibrium: 


0 


cN 


cy  cx 


0 


Constitutive  relation: 
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y 

y 
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0 

N  . 

■V-17 

(27) 


(28) 
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Moment  equilibrium: 


2b 


APPEND  DC  II 


COMPATIBILITY  OF  THE  RESULTANT  LOADS  N  ,  N  .  AND  N 
-  x  ±  y  J. -  xy 

From  Appendix  I,  :  equilibrium  conditions  for  the  x-  and  y-directions 

are  given  by  equations  (27)  as 
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A  constitutive  relation  is  given  by  equation  (8)  as 
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and  equations  (11 )  and  (12)  give 
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(Continued) 


Substitution  of  equations  (11)  and  (12)  into  (6)  and  subsequent  sub 
stitution  of  the  result  into  equations  (27)  yields  equations  of  the 
form 


Liu  +  l2v  -  V  -  Vw  -  Vyz 
+  V  *  Liow  -  '  Vyz 


(32) 


where  the  "L's"  are  linear  operators  defined  as  follows: 
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A66  .  2  +  ^2 6  T~”“  +  A22  2 

ox  u  oxdy  ^  ay 
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Solution  of  the  simultaneous  equations  (32)  gives  for  the  reference- 
surface  displacements 


*  -  (l6l9  -  ¥*)>  w  -  <L6h  -  hty  ~xz  -  (L6L5  -  VbJ 


AV  .  (L^  -  L,^)  w  -  (LlL5  -  1^)  7XZ  -  (L^  .  L2I  ) 


y  yyz 


where 


A  -  I^Lg  -  L2 


APPENDIX  III 


REFERENCE -SURFACE  STRAIN  COMPATIBILITY 

Elimination  of  the  reference  surface-displacements  between  the  strain- 
displacement  relations  given  by  equations  (11 )  results  in  the  following 
strain  compatibility  expression: 
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(35) 


Consideration  of  the  equilibrium  conditions  for  the  x-  and  y-directions 
given  by  equation  (18)  allows  the  definition  of  the  Airy  function 
<t(x,y)  such  that 
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With  the  constitutive  relation  given  by  equation  (13)  and  the  relations 
(36),  the  compatibility  equation  (35)  becomes 
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APPENDIX  IV 


BUCKLING  LOAD  CALCULATION 
DIRECT  VARIATIONAL  METHOD 

Substitution  of  the  assumed  trigonometric  functions  given  by  equations 
(22),  (23),  and  (24)  into  the  functional  (17)  and  use  of  the  following 
integrations 

L 

f  .  2  mjtx  ,  L 

J  Sin  —  d*  =  2 
0 
L 

f  .  mnx  nwx  .  _ 

I  sin  — =—  cos  -r—  dx  =  0 

0 

L 

f  2  mnx  .  L 
J  cos  -j-  dx  =  ? 

0 

/-“'■■I1-  h 

0 

L 

f  nwx  .  _ 

J  cos  — —  dx  =  0 

0 

as  well  as  corresponding  integrations  with  respect  to  y ,  yields 
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The  variation  of  U"  +  V"  with  respect  to  the  free  parameters  gives 


■  | 

6«(u" +  v")  -  f  [  4  ;f  4  4  if ) ♦  <6  if  f  if 2 
+  24iffif)2  -»x  (f  1  » 

o  ' 

+  -  4  iff  ■  2a66  (f  )(f  f  -  4  (f  f  (f )  k 

3  2  2 

+  -4iF-=4ifr,f)-4if;  if:  r, 

+  -4Mff  ^(ff|*i*|-««(f)(f^ 

+  dll  Iff  +  *21  Iff  j  *3  ]  -  o 


22  lb 


(39a) 


Lb 


5  (U"  +  V")  =  r 


Mr) -^fi(r  (rr-4--^1 


12  \b 


nit  \  /mitv 
L  /I 


W 


G  h  +  d., 
x  11 


r 

1— f 

+  d*  1  s«  f 

T  + 

U  I 

66  \  b  1 

X 

*  /  nit  \  ,  mit  >  ,  /  mit  '  /  nit  \ 

Mr\r  +  MrAr) 


-  *12 (r)  M  Mr) 


■V 


j  I  ' 

-Mr; 


N„ 


=  0 


(39b) 


Lb 


B  (U-  +  V'0  =  r 


-4  Iff- 


*  /  ns  \  /  mjt  \  t  mit  \ ,  nit  \ 

Mr  lr)  +  Mr.lr, 

‘66  If) 


V 


*  ,  mit  ,  nn  nit 

2d6«  \r  )r  -Mr  r; 

2  2 
r  .  ■  n«  ,  ,  *  ;  m  ■ 

°vh+<fe,  r  +  d66  r 


-  ^(f  )ri + 


*>♦ 


,,  I  221 

-  *21  [r  * 


N  , 


=  0 


W 


y 

(39c) 


5  (U"  +  V")  =  o 


(39d) 


6  (U"  +  V")  = 


Lb 

r 


(f)vifh 


=  0 


(39e) 


6_  (U"  +  V") 

1 


Lb 

r 


-(fK  +  ifh 


=  o 


(39f) 


6  (U"  +  V")  =  0 

v2 


(398) 


32 


6^(11"  +  V")  =  r 


^(ff +  42i(rf|w“1iilr) 
d21  (r)  ry  '  (f1)  "2  "  *12*1  *  *11*3 


(S9h) 


Lb 


\&”  +  V")  =  f|  -I+d12 


lrf-Mrf|M-Mr)  r. 

"  ^2  (b2)  ry  '[T )  V1  '  a22Nl  "  *12N3 


0 

(391) 


6n2(U-  +  V") 


Lb 

r 


=  0 
(39J) 


Through  the  use  of  equations  (39e),  (39f)>  (39h),  and  (391 )  to  eliminate 
N3,  N^,  N2,  and  v^,  the  system  of  equations  (39)  is  reduced  quite  readily 
to  a  set  of  linear,  homogeneous  equations  in  terms  of  W,  r^,  and  Ng* 
A  non-trivial  solution  of  the  set  of  f  ur  equations  exists  only  if  the 
determinant  of  the  coefficients  vanishes;  that  is, 
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where  the  various  nondimens ional  parameters  are  defined  as  follows: 
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APPENDIX  V 


VIBRATION  PROBLEM 

The  governing  equations  for  the  free  lateral  vibration  of  the  circular 
cylindrical  shell  can  be  derived  using  the  following  form  of  Hamilton's 
variational  principle: 

s  h 

I  (U"  -  T ")  dt  -  0  (12) 

-tl 

In  equation  (12),  U"  represents  the  Reissner  form  of  the  strain 
energy  and  T"  represents  the  kinetic  energy.  In  this  case,  both 
U"  and  T"  are  functions  of  time.  With  the  assumption  of  harmonic 
motion,  the  time  dependence  of  equation  (12)  is  eliminated  and  the 
variational  principle  is  stated  as 

5(U"  -  T")  =  0  (13) 


For  the  fiber-reinforced  shell,  U"  '  *  given  in  equation  (IT ) .  When 
the  radial  (lateral)  motion  of  the  shell  is  assumed  to  predominate, 

T"  is  given  by 

2 
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iji  /  /  _  _  __ 

2 

where  u  is  the  natural  frequency  and  p  is  the  mass  density  of  the 
shell  per  unit  reference  surface  area.  The  functional  for  the  vari¬ 
ational  procedure  is  expressed  as 
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Variation  of  the  functional  given  by  equation  (45)  with  respect  to 
Nx  ,  Ny  ,  ,  u  ,  v  ,  w  ,  7^^  ,  ~yz  yields  the  Euler  equations  and 

associated  boundary  conditions  for  the  vibration  problem.  If  Nx  is 
set  equal  to  zero,  the  boundary  conditions  and  the  Euler  equations, 
except  for  the  lateral  equilibrium  equation,  are  identical  to  those 
established  in  Appendix  I  for  the  buckling  problem.  The  lateral 
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equilibrium  equation  for  the  vibration  problem  becomes 
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With  the  use  of  the  assumed  generalized  displacements  and  forces  given 
by  equations  (22),  (23),  and  (2b)  with  set  equal  to  zero,  the 

direct  variational  procedure  yields  a  frequency  equation  in  determinant 
form  as 
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APPENDIX  VI 

DEFINITIONS  OF  MATRICES 
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[a]  = 

[A]”1 

(51) 

[d']  = 

[  a]  [D] 

(52) 

[d]  = 

[D ]  [ a ]  =  [d'  ]T 

(53) 

[d*]  = 

[D*]-[D][a][D]  =  [D4]-[D][d7] 

(54) 
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The  effects  of  transverse  shear  deformations  upon  the  classical  buckling  load  of 
fiber-reinforced  cylindrical  shells  under  uniform  axial  compression  have  been 
analyzed.  A  stability  determinant  which  upon  evaluation  yields  an  expression  for 
the  buckling  load  has  been  developed  using  a  modified  form  of  the  Reissner 
variational  principle.  The  buckling  loads  predicted  by  the  stability  determinant, 
with  transverse  shear  effects  neglected,  have  been  calculated  and  shown  to  agree 
with  previously  published  results  for  boron-epoxy  and  glass-epoxy  cylinders. 
Inclusion  of  the  transverse  shear  effects  in  the  two  cases  investigated  shows 
little  reduction  in  the  classical  buckling  loads.  For  general  application,  charts 
are  presented  to  give  stability  criteria  for  fiber-reinforced  composites  as  a 
function  of  the  geometric  and  mechanical  properties.  Effective  transverse  shear 
rigidities,  established  by  experiment,  must  be  developed  in  order  to  realistically 
estimate  the  reduction  in  buckling  load  from  the  charts  presented.  For  the 
determination  of  natural  frequencies  of  fiber-reinforced  shells,  with  transverse 
shear  effects  included,  a  frequency  equation  in  determinant  form  is  presented. 
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